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psj ■ Abstract 



A powerful method to solve nonlinear first-order ordinary differential equations, 
which is based on a geometrical understanding of the corresponding dynamics of 
the so-called Lie systems, is developed. This method enables us not only to solve 
some of these equations, but also gives geometrical explanations for some, already 
known, ad hoc methods of dealing with such problems. 
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1 Introduction 

^ ■ Systems of nonautonomous first-order ordinary differential equations appear often in 

. . ! Mathematics, Physics, Chemistry and Engineering and therefore methods to solve them 

. ^ \ and analyse their properties are specially interesting because they allow us to understand 

^ \ many important problems in these various fields. 

■ A special kind of these systems, the so-called Lie systems (or Lie-Scheffers systems) 

[l]-[5], has recently been analysed in many papers [6]-[Tl]. An important property of 
such systems is that they admit a certain set of time-dependent changes of variables 
which transforms each Lie system into a new one [TO]. Such changes can be used to 
transform a given Lie system into an easily integrable one, i.e. into a Lie system related 
to a solvable Lie algebra of vector fields. In all these cases, we can obtain constants of 
the motion, integrability conditions or even solutions [T^. 

This transformation property is not only valid for Lie systems but it still holds for a 
more general set of systems of differential equations. In order to generalise this property 
we develop the concept of a quasi-Lie system, i.e. a system of nonautonomous differential 
equations accompanied with a flow of diffeomorphisms transforming it into a Lie system. 
We must stress that quasi-Lie systems are not equivalent to Lie systems in the trivial 
sense: the transformations we use are time-dependent and we might get time-dependent 
superposition rules for them, meanwhile. Lie systems are related to superposition rules 
which are time-independent. Furthermore we define quasi-Lie schemes. A scheme is a 
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procedure that, in some cases, enables us to transform a given nonautonomous system of 
differential equations into a Lie system. This scheme, when applied to Lie systems, gets 
back to the well-known transformations of Lie systems, see [7]. However, many other 
important systems of first-order differential equations can be studied by means of this 
new method. In this sense the method of quasi-Lie schemes may be viewed as an abstract 
unification of many ad hoc developed integration methods existing in the literature. 

The aim of this paper is to introduce and to analyse the properties of quasi-Lie 
schemes and to illustrate the theory by developing some interesting examples. More 
specifically, in order to the paper be self-contained. Section 2 is devoted to review the 
theory of time-dependent vector fields. Then we summarise in Section 3 some prop- 
erties of Lie systems and in Section 4 we study the quasi-Lie scheme concept which 
may be used to relate certain sort of systems of differential equations to Lie systems. 
Such quasi-Lie schemes are applied in Section 5 to analyse some interesting systems of 
differential equations. Some well-known differential equations are analysed: dissipative 
Milne-Pinney equations, nonlinear oscillator and Emden equations in order to recover 
from an unified point of view some of their properties. The fields of applications of these 
differential equations is very broad, i.e. just the Emden equation appears in Mathemat- 
ical Physics [15], Theoretical Physics [TB], Astronomy [T7], Astrophysics [IB] and the 
review by Wong [TH] contains about 140 references. We also study a time-dependent 
dissipative Mathews-Lakshmanan oscillator and we provide a, as far as we know, new 
time-dependent constant of the motion. Finally, in Section 6, we sum up the conclusions 
of our paper and give an outlook of possible problems to be studied by means of the 
methods here developed. 



2 Generalised flows and time-dependent vector flelds 

A nonautonomous system of first-order ordinary differential equations in a manifold 
is represented by a time-dependent vector field X = X{t,x) on such a manifold. On 
a non-compact manifold, the vector field Xt{x) = X{t,x), for a fixed t, is generally not 
defined globally, but it is well defined on a neighbourhood of each point xq E N for 
sufficiently small t. It is convenient to add the time t to the manifold and to consider 
the autonomisation of our system, i.e. the vector field 

X{t,x) = ^ + 

defined on a neighbourhood A^"^ of {0} x in M x A^. The vector field Xt is then defined 
on the open set of A^, 

= {x,eN\{t,x^)eN''}, 

for all t G M. If = N for all t G M, we speak about a global time- dependent vector 
field. The system of differential equations associated with the time- dependent vector 
field X{t,x) is written in local coordinates 

—— = X^(t,x), i = 1, . . . ,n = dim N, 
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where X{t,x) = ^"^^ is defined locally on the manifold for sufficiently 

small t. 

A solution of this system is represented by a curve s t— > 7(3) in N (integral curve) 
whose tangent vector 7 at t, so at the point 7(t) of the manifold, equals X(t,'j(t)). In 
other words, 

7(t)=X(t,7(t)). (1) 

It is well-known that, at least for smooth X which we work with, for each xq there is 
a unique maximal solution 7^°(t) of system ([T]) with the initial value Xq, i.e. satisfying 
7x (0) ~ ^0- This solution is defined at least for t's from a neighbourhood of 0. In case 
7^°(t) is defined for all t G M, we speak about a global time solution. 

The collection of all maximal solutions of the system ([1]) gives rise to a (local) gen- 
eralised flow on A^. By a generalised flow g on N we understand a smooth time- 
dependent family gt of local diffeomorphisms on A^, gt{x) = g(t,x), such that go = id^. 
More precisely, g is a. smooth map from a neighbourhood of {0} x N mM.x N into N, 
such that gt maps diffeomorphically the open submanifold = {xq & N \ {t, xq) G N^} 
onto its image, and go = idjv- Again, for each xq G A^ there is a neighbourhood 
of Xq in A^ and e > such that gt is defined on U^o for t G (— e, e) and maps Ux:^ 
diffeomorphically onto gtiU^^^)- 

If A^f = A^ for all t G M, we speak about a global generalised flow. In this case 
g -.t eM.^ gt & Diff(A^) may be viewed as a smooth curve in the diffeomorphism group 
Diff(A^) with ^0 = idiv- 

Here is also convenient to autonomise the generalised flow g extending it to a single 
local diffeomorphism 

g{t,x) = {t,g{t,x)) (2) 

defined on the neighbourhood A^^ of {0} x A^ in M x A^. The generalised flow g^ induced 
by the time-dependent vector field X is defined by 

<7^(t,xo)=7x"W. (3) 
Note that, for g = g^ , equation ([3]) can be rewritten in the form: 

Xt = X{t,x)=gtog-\ (4) 

In the above formula, we understood Xt and gt as maps from A^ into TN, where gt{x) 
is the vector tangent to the curve s ^— > g{s,x) at g{t,x). Of course, the composition 
9t ° 9t^) called sometimes the right-logarithmic derivative oi t ^ gt, is defined only for 
those points Xq G A^ for which it makes sense. But it is always the case locally for 
sufficiently small t. 

Let us observe that equation (jl]) defines, in fact, a one-to-one correspondence be- 
tween generalised flows and time-dependent vector fields modulo the observation that 
the domains of gt o g^^ and Xt need not coincide. In any case, however, gt o g^^ and 
Xt coincide in a neighbourhood of any point for sufficiently small t. One can simply say 
that the germs of X and gt o g^^ coincide, where the germ in our context is understood 
as the class of corresponding objects that coincide on a neighbourhood of {0} x A^ in 
M X A^. 
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Indeed, for a given g, the corresponding time-dependent vector field is defined by (jlj). 
Conversely, for a given X, the equation (jlj) determines the germ of the generalised flow 
g{t, x) uniquely, as for each x = Xq and for small t equation (jl]) implies that t \^ g(t, xq) 
is the solution of the system defined by X with the initial value xq. In this way we get 
the following 

Theorem 1. Equation ^ defines a one-to-one correspondence between the germs of 
generalised flows and the germs of time- dependent vector fields on N. For compact N, 
this correspondence reduces to a one-to-one correspondence between global time- dependent 
vector fields and global generalised flows. 

Any two generalised flows g and h can be composed: by definition {g o h)t = gt ° ht, 
where, as usual, we view gt oht as a. local diffeomorphism defined for points for which the 
composition makes sense. It is important that in a neighbourhood of any point it really 
makes sense for sufficiently small t. As generalised flows correspond to time-dependent 
vector flelds, this gives rise to an action of a generalised flow h on a time-dependent 
vector fleld X, giving rise to h^X deflned by the equation 

= hog^. (5) 

To obtain a more explicit form of this action, let us observe that 

{h^X)t = o (h o g^);^ = {ht o g^ + Dhtm) o (g^)-- o h^\ 

and therefore 

{h^X)t = ht o h-^ + Dht [gf o {g^)-^) o h-\ 

i.e. 

{h^X)t = htoh^' + {ht).{Xt), (6) 

where (/i*)* is the standard action of diffeomorphisms on vector flelds. In a slightly 
different form, this can be written as an action of time-dependent vector flelds on time- 
dependent vector flelds: 

{9lX)t = Yt + {gJUXt)- (7) 
For global time- dependent vector flelds on compact manifolds the latter deflnes a group 
structure in global time-dependent vector flelds. This is an inflnite-dimensional analog 
of a group structure on paths in a flnite-dimensional Lie algebra which has been used 
as a source for a nice construction of the corresponding Lie group in [20] . Since every 
generalised flow has the inverse, {g~^)t = {gt)~^, so generalised flows, or better to say, the 
corresponding germs, form a group and the formula ([7]) allows us to compute the time- 
dependent vector fleld (right-logarithmic derivative) Xj~^ associated with the inverse. It 
is the time-dependent vector fleld 

Xi' = -{gf):\Xt) . (8) 

For time- independent vector flelds Xt = Xq for all t we have {g^)^X = X and also we 
get the well-known formula 

x-^ = -X. 

Note that, by deflnition, the integral curves of h^X are of the form ht{'j{t)), where 7(t) 
are integral curves of X. We can summarise our observation as follows. 
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Theorem 2. The equation ^ defines a natural action of generalised flows on time- 
dependent vector fields. This action is a group action in the sense that 

(g o h^X = g^ih^X) . 

The integral curves of h^X are of the form ht{'y{t)), for'y{t) being an arbitrary integral 
curve for X . 

The above action of generahsed flows on time-dependent vector flelds can be also deflned 
in an elegant way by means of the corresponding autonomisations. It is namely easy to 
check the following. 

Theorem 3. For any generalised flow h and any time- dependent vector field X on a 
manifold N, the standard action h^,X of the diffeomorphism h, being the autonomisation 
of h, on the vector field X , being the autonomisation of X, is the autonomisation of the 
time- dependent vector field h^X: 

KX = Kpc. 

3 Lie systems and superposition rules 

The conditions for the system determined by a time-dependent vector fleld X{t, x) on 
a manifold ensuring that it admits a superposition rule, i.e. that there exists, at least 
locally for an open U C iV™ x M", a map $ : f/ ^ N, x = . . . , X(rn)', ki, . . . , kn), 

such that its general solution can be written as 

x{t) = $(x(i)(t),...,X(„)(t);/ci,...,/c„), 

where | a = 1, . . . ,m} is any family of particular solutions 'in general position' 

and k = {ki, . . . , kn) is a set of n arbitrary constants, were studied by S. Lie [T]. Let us 
stress that the superposition function $ is time-independent. 

The necessary and sufficient conditions say that the associated time- dependent vector 
field X can be written as a linear combination 

r 

Xt = Y,ba{t)X^o.) , (9) 

a=l 

such that the vector fields | a = 1, . . . ,r} generate a finite-dimensional real Lie 

algebra, the so-called Vessiot-Guldberg Lie algebra. The latter means that there exist 
real numbers Ca/s-y, such that 

r 

[X(„) , X(^)] = ^ Ca/3^X(^) , a, P = 1, . . . ,r . 

7=1 

Linear systems are particular instances of Lie systems associated with a Vessiot-Guldberg 
Lie algebra, isomorphic to the Lie algebra Qi{n, M), and m = n in the homogeneous case. 
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or the corresponding affine algebra and m = n + 1 in the inhomogeneous one. The 
Riccati equation is another example for which Xf = 6o(^)^o + bi{t)Yi + &2(^)^2 with 

ox ox ox 

closing on a Vessiot-Guldberg Lie algebra isomorphic to the Lie algebra s[(2, M), see [71(9]. 
There is an action <l>Ric : SL{2, M) x 1 ^ M of the Lie group SL{2, i?) on 1 = M U cx) ^ 
such that the fundamental vector fields of this action are linear combinations with real 
coefficients of the vector fields Yq, Yi and 12- 

Another relevant example of a Lie system is given by a time-dependent right-invariant 
vector field in a Lie group G. If {ai, . . . , a^} is a basis of T^G and are the corre- 
sponding right-invariant vector fields, X^{g) = {Rg)^a.a, then the time-dependent right- 
invariant vector field 

r 

x^ = -Y,K{t)Xl 

a=l 

defines a Lie system in G whose integral curves are solutions of the system g = 
— '^a=iba{t)X^{g). Applying to both sides, we see that g(t) satisfies 

r 

{Rg-i(t))*g{t) =-Y, ba{t)aa G TeG . (10) 

a=l 

Right-invariance means that it is enough to know one solution, for instance the one 
starting from the neutral element e, to know all the solutions of the equation with any 
initial condition, i.e. we obtain the solution g'{t) with the initial condition g'{0) = go as 
Rg{o)g{t)- A generalisation of the method used by Wei and Norman for linear systems [21] 
is very useful in solving such equations and furthermore there exist reduction techniques 
that can also be used [TT]. Finally, as right-invariant vector fields X^ project onto the 
fundamental vector fields in each homogeneous space of G, the solution of ffTUl) allows 
us to find the general solution for the corresponding Lie system in each homogeneous 
space. Conversely, the knowledge of particular solutions of the associated system in a 
homogeneous space gives us a method for reducing the problem to the corresponding 
isotropy group pAj. This equation is also important because any Lie system described 
by a time-dependent vector field on a manifold A^, like (12]), where the vector fields are 
complete and satisfy the same commutation relations as the basis {ai, . . . , a.r}, determines 
an action $LieG '■ GxN N such that the vector field Xa is the fundamental vector field 
corresponding to a.a, and moreover, the integral curves for the time-dependent vector field 
are obtained from the solutions of equation flTUl) . More explicitly, the general solution 
of the given Lie system is x{t) = ^hieG{git)yXo)j where Xq is the initial condition of the 
solution and g(t) is the solution of equation flTU]) with g{0) = e. 

The search for the number m of solutions and the superposition function $ has 
recently been studied from a geometric perspective [12]. Essentially, we should consider 
'diagonal prolongations' to A^"^"*"^, X{x(o), . . . , i), of the time- dependent vector field 
X{t,x) = ^^=iX^{t,x) d/dx^, given by 

m 

X(a;(o), . . . = ^Xa(a;(a),t) , t G M, 

a=0 
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where Xa{x(^a),t) = Sr=i ^^^l-^W' ^) ^/^-^(a)' ^^^^ ^^^^ extended system admits n 
independent constants of the motion, which define in an imphcit way the superposition 
function. 

4 Quasi-Lie systems and schemes 

Definition 1. By a quasi-Lie system we understand a pair {X,g) consisting of a time- 
dependent vector field X on a manifold N (the system) and a generalised fiow g on N 
(the control) such that QifX is a Lie system. 

Since for the Lie system g-^^X we are able to produce the general solution out of a number 
of known particular solutions, the knowledge of the control makes a similar procedure for 
our initial system possible. Indeed, let $ = $(xi, . . . , Xm, ki, . . . , kn) be a superposition 
function for the Lie system g^^X, so that, knowing m solutions X{i), . . . , x^m) of Qi^X, we 
can derive the general solution of the form 

X(o) = $(X(i), . . . , fcl, . . . , kn) . 

If we now know m independent solutions, . . . , of X, then, according to Theorem 
El Xait) = gt{xa(t)) are solutions of g^^X, producing a general solution of g^X in the 
form $(a;(i), . . . ^i, • • • , ^n)- It is now clear that 

X{0)it) = gt^ O ^{gt{x(^i){t)), . . . ,gt{X(^rn){t)), fci, . . . , kn) (11) 

is a general solution of X. In this way we have obtained a time- dependent superposition 
rule for the system X. We can summarise the above considerations as follows. 

Theorem 4. Any quasi-Lie system {X,g) admits a time- dependent superposition rule 
of the form / flT]) . where ^ is a superposition function for the Lie system gifX . 

Of course, the above time-dependent superposition rule is practically meaningful in find- 
ing the general solution of a system X only if the generalised fiow g is explicitly known. 
An alternative abstract definition of a quasi-Lie system as a time-dependent vector field 
X for which there exists a generalised fiow g such that g^^X is a Lie system does not 
have much sense, as every X would be a quasi-Lie system in this context. For instance, 
given a time-dependent vector field X, the pair (X, {g^)~^) is a quasi-Lie system because 
{g'^)^^ ° gf = idAT, thus {g^)^^X = 0, which is a Lie system trivially. On the other 
hand, finding (f?^)"^ is nothing but solving our system X completely, so we just reduce 
to our original problem. In practice, it is therefore crucial that the control g comes 
from a system which can be integrated effectively. There are, however, many cases when 
our procedure works well and provides a geometrical interpretation of many, originally 
developed ad hoc, methods of integration. Consider, for instance, the following scheme 
that can lead to 'nice' quasi-Lie systems. 

Take a finite-dimensional real vector space V of vector fields on N and consider the 
family of all time-dependent vector fields X on X such that Xt belongs to V on its 
domain, i.e. Xt e V^n^- We will say that these are time-dependent vector fields taking 
values in V. The time-dependent vector fields taking values in V depend on a finite family 
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of control functions. For example, take a basis {Xi, . . . , Xr} of V and consider a general 
time-dependent system with values in V determined by 6 = b{t) = . . . , br(t)) as 

{X'), = b,{t)X, . 

On the other hand, the nonautonomous systems of differential equations associated with 
X e V\]^x are not Lie systems in general, if V is not a Lie algebra itself. If we have 
additionally a finitely parametrised family of local diffeomorphism, say g = g{ai, . . . ,ak), 
then any curve a = a{t) = (ai(t), . . . , akit)) in the control parameters, defined for small 
t, gives rise to a generalised flow g^ = g{a{t)). Let us assume additionally that there 
is a Lie algebra Vq of vector fields contained in V. We can look for control functions 
a{t) such that for certain b{t) we get that g'^X^ has values in Vq for each time t. Let us 
denote this as 

glX' e Vo. (12) 

This choice of control functions makes {X^, g"-) into a quasi-Lie system, so we get time- 
dependent superposition rules for the corresponding systems X^. 

Let us observe that in the case when all the generalised flows g"' preserve V, i.e. for 
each time-dependent vector field X'^ e V also g'^X'^ G V, the inclusion (fT2!) becomes a 
differential equation for the control functions a{t) in terms of the functions b{t). This 
situation is not so rare as it may seem to be at the first sight. Suppose, for instance, 
that we find a Lie algebra W <Z V such that [PV, V^] C and that the time-dependent 
vector fields with values in W can be effectively integrated to generalised flows. In this 
case, any time- dependent vector field with values in W gives rise to a generalised 
flow g"' which, in view of transformation rule ([7]), preserves the set of time-dependent 
vector fields with values in V . For each b = b{t) the inclusion f|T2|) becomes therefore 
a differential equation for the control function a = a{t) which often can be effectively 
solved. 

Definition 2. Let W, V be finite-dimensional real vector spaces of vector fields on a 
manifold A^. We say that they form a quasi-Lie scheme S{W, V) if the following are 
satisfied 

1. W is a. vector subspace of V. 

2. W is a. Lie algebra of vector fields, i.e. [W, W] C W. 

3. W normalises V, i.e. [W, V] C V. 

If is a Lie algebra of vector fields V, we call the quasi-Lie scheme S(y, V) simply a 
Lie scheme S(y). 

Remark 1. There is the largest Lie subalgebra we can use as W - the normalizer of V 
in V. Sometimes, however, it is useful to consider smaller Lie subalgebras W. 

We say that a time-dependent vector field X is in a quasi-Lie scheme S(W, V), and 
write X G S(W,V), if X belongs to V on its domain, i.e. Xt G V^n^- Note that, by 
definition, the set of all time-dependent vector fields belonging to S{W, V) depends only 
on V and the choice of W is irrelevant. 

Now, given a quasi-Lie scheme S{W, V) which we will call sometimes simply a scheme, 
we may consider the group Q(W) of generalised flows associated with W. 
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Definition 3. We call the group of the scheme S{W, V) the group Q{W) of generalised 
flows corresponding to the time-dependent vector fields with values in W . 

Proposition 1. (Main property of a scheme) Given a scheme S{W,V), a time- 
dependent vector field X G S(W,V), and a generahsed flow g G GiW), we get that 
g^XeS{W,V). 

The proof is obvious and follows from the fact that if y is a generalised flow with 
values in W and X takes values in V, then, according to the formula ([7]), g^X takes 
values in V as well, as [W, V] (ZV and V is finite-dimensional. 

From the last definition we can state the definition of quasi-Lie system with respect 
to a scheme. 

Definition 4. Given a quasi-Lie scheme S{W, V) and a time-dependent vector field 
X G S(W, V), we say that X is a quasi-Lie system with respect to S{W, V) if there exists 
a generalised flow g G Q{W) and a Lie algebra of vector fields Vq G V such that 

g^X G SiVo). 

We emphasise that if X is a quasi-Lie system with respect to the scheme S{W, V), it 
automatically admits a time- dependent superposition rule in the form given by (fTTl) . In 
the next Section, we apply our theory and illustrate these concepts with examples. 

5 Applications of quasi-Lie schemes 

The above mentioned properties of quasi-Lie schemes and quasi-Lie systems can be 
used to investigate some previously studied systems of differential equations [22]- [33] 
systematically from this new perspective. 

In this Section, we apply quasi-Lie schemes to study dissipative Milne-Pinney equa- 
tions, nonlinear oscillators and Emden differential equations. More precisely, we first ap- 
ply our theory to dissipative Milne-Pinney equations. These systems cannot be treated 
with the theory of Lie systems directly but one can use a quasi-Lie scheme to transform 
them into Milne-Pinney equations. These latter equations have been proved to be SODE 
Lie systems recently [6], [13] and this fact enables us to get a time-dependent superpo- 
sition rules for dissipative Milne-Pinney equations by means of the superposition rule 
found for non-dissipative ones. 

Next, we analyse nonlinear oscillators. Perelomov studied some of these systems 
in order to relate them to other important systems [2S]. The cases investigated by 
Perelomov were selected and obtained by means of ad hoc methods. Here we approach 
some instances treated in [28] in order to explain Perelomov's work from the point of 
view of the theory of quasi-Lie schemes. As a result, we show how our theory provides 
time-dependent constants of the motion and clarify some points about this work. 

We also analyse Emden equations. In this case, we use a quasi-Lie scheme to obtain 
constants of the motion for Emden equations whose time- dependent coefficients satisfy 
certain conditions. Notwithstanding, this is just a small instance of what can be made 
by means of our scheme. 



9 



Finally, we analyse certain kind of dissipative Mathews-Lakshmanan oscillators 
121 ES]- Some kinds of these nonlinear oscillators have been recently investigated 
from the point of view of Classical and Quantum Mechanics [331 [311 ES]- Here we just 
perform a simple application of the theory of quasi-Lie schemes to relate different types 
of dissipative Mathews-Lakschmanan oscillators to the usual ones. 

Let us provide a quasi-Lie scheme to deal with some of the systems investigated in 
following subsections. Recall that we need to find vector spaces W and V of vector 
fields satisfying the three conditions stated in definition [2J Consider the vector space V 
spanned by the linear combinations of the vector fields 



d 



X- 



d 



X 



X, 



d 



d 



OV OV OX OV 

on TM and take the vector subspace W <Z V generated by 



X, 



dx 



(13) 



n = X4 = V—, 



d_ 

dv' 



Xx = x-^, 

OV 



d 

F3 = X5 = X — . 

ox 



Therefore, is a solvable Lie algebra of vector fields, 

[n,i2] = -r2, [1^1,^3] = 0, [y2,i3] = -y2, 

and taking into account that 



[YuX^] 



-X2, 

X5 — X4, 



[n,X3] 
[>3,X2] 



X3, 

nX2, 



[>^2,X2] 

[Y,.X,] 



0, 



-X3 



we see that V is invariant under the action of W , i.e. V] C V . In this way we get 
the quasi-Lie scheme S{W,V). We stress that the vector space V is not a Lie algebra 
because the commutator [X2, X3] is not in V . Moreover, there is no Lie algebra of vector 
fields V and thus V cannot be related to a Lie scheme. 

The key tool provided by the scheme S{W, V) is the infinite-dimensional group Q{W) 
of generalised flows for the time-dependent vector fields with values in W , i.e. ai{t)Yi + 
Oi2{i)Y2+a^{t)Y^. The integration of such t-dependent vector fields lead to the description 
of the time-dependent changes of variables associated with Q{W), i.e. 



g{W) = {g{a{t),m.l{t)) 



V = a{t)v' + I3{t)x' 



a(t),7W > 0, a(0) = 7(0) = 1, /3(0) = 



5.1 Dissipative Milne— Pinney equations 

In this Section, we study the so-called dissipative Milne-Pinney equations. We 
show that the first-order ordinary differential equations associated with these second- 
order ones in the usual way, i.e. by considering velocities as new variables, are not 
Lie systems. However, the theory of quasi-Lie schemes can be used to deal with such 
first-order systems. Here we provide a scheme which enables us to transform a certain 
kind of dissipative Milne-Pinney equations, considered as first-order systems, into some 
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first-order Milne-Pinney equations already studied by means of tlie tfieory of Lie systems 
[0]. As a result we get a time- dependent superposition rule for some of these dissipative 
Milne-Pinney equations. 

Let us state the problem under study. Consider the family of dissipative Milne- 
Pinney equations of the form 



We are mainly interested in the case c{t) 7^ 0, so we assume that c(t) has a constant sign 
for the set of values of t we are considering. 

Usually, we associate with such a second-order differential equation a system of first- 
order differential equations by introducing a new variable v and relating the differential 
equation f|T^ to the system of first-order differential equations 



In order to verify that the quasi-Lie scheme S(W,V), for the case n = —3, can be 
used to handle the latter system, that is, X G S{W,V), we have to ensure that the 
time-dependent vector field 



whose integral curves are solutions for the system (ITS!) . is such that Xt & V for every t 
in an open interval in M. In this way, in view of (1131) . we observe that 



and thus X G S{W, V). Moreover V" = {Xi, . . . , X4) is not a Lie algebra of vector fields 
because [X3, X2] ^ V". Also, there is no finite-dimensional real Lie algebra V containing 
V". Thus, system (fT5|) is not a Lie system but we can use the quasi-Lie scheme S(W, V) 
to investigate it. 

Let us consider the infinite-dimensional subgroup of Q{W) given by its time-dependent 
changes of variables with 7(t) = 1. According to the general theory of quasi-Lie schemes, 
these time-dependent changes of variables enable us to transform system f|T5|) into a new 
one again describing the integral curves for a time-dependent vector field X' G S{W, V), 
that is. 



X = a{t)x + h{t)x + c(t)— . 



(14) 




(15) 




Xt = a{t)X^ + 6(t)Xi + c{t)X2 + X3, 



X[ = a'(t)X4 + b'it)Xi + c'{t)X2 + d'{t)Xs + e\t)X^ . 



(16) 



The new coefficients are 



a'{t) 



ait) - m 



a{t) 




+ ait) 



a{ty 

m p'it) m 



< 



a(t) a(t) a(t)^ 



d'{t) = a{t), 
.e'{t)=m- 
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The integral curves for the time-dependent vector field f|T6|) are solutions of the system 



^ = f3[t)x' + a{t)v', 

dt \a{t) ^ ^a{t) a{t) a{t) ' ^ ' 

( , ^ n, ^ <^{t)\ , c(t) 1 
+ a{t) - (3{t) - -7-( N + 



a{t) J a{t)x'^' 

As it was said in Section HJ we use schemes to transform the corresponding systems 
of first-order differential equations into Lie ones. So, in this case, we must find a Lie 
algebra of vector fields Vq (Z V and a generalised fiow g G Q(W) such that g^X G 5'(Vo). 
This leads to a system of ordinary differential equations for the functions a{t), (3{t) and 
some integrability conditions about the initial functions a{t),b{t) and c{t) for such a 
time-dependent change of variables to exist. 

In order to find a proper Lie algebra of vector fields Vq C V, note that Milne-Pinney 
equations studied in [6] are Lie systems in the family of differential equations defined by 
systems (1151) and therefore it is natural to look for the conditions needed to transform 
a given system ( fT5i) . described by the time- dependent vector field Xt, into one of these 
first-order Milne-Pinney equations of the form 

-uj{t)x + f{t) — , 

where is a constant, i.e. a system describing the integral curves for a time-dependent 
vector field with values in the Lie algebra of vector fields [6] 

Vo = {X, + kX2,X^,^{X,-X,)). 

As a result, we get that P = 0, a = f and, furthermore, the functions a, a and c must 
satisfy 

ka"^ = c, a — aa = 0, (19) 

so that c must be of constant sign equal to that of k. The second condition is a differential 
equation for a and the first one determines c in terms of a. Therefore, both conditions 
lead to a relation between c and a providing the integrability condition 



c{t) = kexp{2A{t)) , and A{t) = J a{t)dt, (20) 

and showing, in view of f|T7j) . ffTSjl and f|T9|) . that 

a{t) = exp {A{t)) and iu{t) = -h(t) exp (-A(t)) , 
where we choose the constants of integration in order to get a(0) = 1 as required. 
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Summing up the preceding results, under the integrabihty condition fl20l) . the first- 
order Milne-Pinney equation f[T^ can be transformed into the system 

^ ^ = exp{A{t))v', 

^ = 6(t)exp(-A(t))x' + exp(A(t)) A, 
k at x'-^ 

by means of the time-dependent change of variables 
^(exp (A(t)),0,l) = 



X 

v' 



exp{A{t))v 



Note 5.1. The previous change of variables is a particular instance of the so-called 
Liouville transformation JEM/- 

Now, the final Milne-Pinney equation can be rewritten by means of the time 
reparametrisation 

r(t) = j exp {A{t))dt, 

as 



dx' , 

- =v, 

^ =exp(-2A(t))6(t(r))x' + 4- 
dr X 



These systems were analysed in [7] and there it was shown through the theory of Lie 
systems that they admit the constant of the motion 



™ X 2 



/ = (xv' — vx'Y + k ( — 
\x' 

where (x, v) is a solution of the system 

dx 

dr ' 

^ = exp(-2A(t))b(t)x, 
dr 

which can be written as a second-order differential equation 

^ = exp{-2A{t))b{t)x. 

If we invert the time reparametrisation, we obtain the following differential equation 

:^ - a(t)l - 6(t)x = 0, (21) 
which is the linear differential equation associated with the initial Milne-Pinney equation. 
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As it was shown in [6], we can obtain, by means of the theory of Lie systems, the 
following superposition rule 

./O / , X 1/2 

/ V ^ / r -2 , r -2 I ATT 1 I- - - - \o - - \ ' 



\XxV2 - VlX2\ 



(^hxl + hxl ± A/4/1/2 - k{xiV2 - V1V2Y X1X2 



and as the time- dependent transformation performed does not change the variable x, we 
get the time-dependent superposition rule 

V^<^{t) f T -2 . r -2 , /, T T ^ _ _ _ ^ 

X = T—^ ^—7 I2X-1 + hxo ± \ 4:Iil2 7rr:r{xiX2 - X1X2Y X1X2 

\x1X2 - X1X2I \ y a'^[t) 

in terms of a set of solutions of the second-order linear system f l21l) . 

Summing up, the application of our scheme to the family of dissipative Milne-Pinney 
equations 

X = a(t) X + b(t) X + exp ^2 J a{t) dt^ — 
shows that it admits a time-dependent superposition principle: 

1/2 



V2a(t) T 2 T 2 ^ / 

X = : rr hVi + /1Z/2 ± W 4/1/2 TTT^iyiy^ - l/2Z/i)^ I/1Z/2 

\y1y2 - 1/22/1 1 Y V OL\t) 

in terms of two independent solutions ?/i,?/2 for the differential equation 

y-a{t)y-h{t)y = Q. 

So, we have fully detailed a particular application of the theory of quasi-Lie schemes 
to dissipative Milne-Pinney equations. As a result, we provide a time-dependent super- 
position rule for a family of such systems. Another paper dealing with such an approach 
to dissipative Milne-Pinney equations and explaining some of their properties can be 
found in [37] . 



5.2 Nonlinear oscillators 

As a second application of our theory, we use quasi-Lie schemes to deal with a 
certain kind of nonlinear oscillators. The main objective of this Section is to explain 
some properties of some time-dependent nonlinear oscillators studied by Perelomov in 
[28] . We also furnish with a, as far as we know, new constant of the motion for these 
systems. 

Consider the subset of the family of nonlinear oscillators investigated in [28] : 

X = h{t)x + c(t)a;", 7^ 0, 1 . 

The cases n = 0, 1, are omitted because they can be handled with the usual theory of Lie 
systems. Like in the above Section, we link the above second-order ordinary differential 
equation to the first-order system 

X = V, , ^ 

(22) 

V = h{t)x + c(t)x". ^ ' 
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Now, we have to go over whether the solutions of system fl22|) are integral curves for 
a time- dependent vector field X G S{W, V). In order to check out this, we realise that 
the system describes the integral curves for the time-dependent vector field 

ox ov 

which can be written as 

Xt = h{t)X^ + c{t)X2 + X^. (23) 

Note also that [X2,X3] ^ V and not only V" = {Xi, X2, X3) is not a Lie algebra of 
vector fields but also there is no finite-dimensional Lie algebra V including V". Thus, 
X cannot be considered as a Lie system and we conclude that the first-order nonlinear 
oscillator fl22l) describing integral curves of the time-dependent vector field fl23|) (which 
is not a Lie system) can be described by means of the quasi-Lie scheme S{W, V). 

Let us restrict ourselves to analyse those time-dependent changes of variables asso- 
ciated with the generalised flows of Q{W) with P{t) = 7(t) and a(t) = l/7(t) and apply 
these transformations to the system fl22l) . The main theorem of the theory of quasi-Lie 
systems tells us that 

gia{t),P{t),^it)UXeSiW,V). 
Indeed, these time-dependent transformations lead to the systems 




(24) 

('iHtm - 'lithit)}^' + <:(tyr\t)x'", 

which are related to the second-order differential equations 

7^(t)x' = -27(t)7(t)i' + il'mt) - i{t)^{t))x' + c{t)^-^\t)x'- . 

But the theory of quasi-Lie schemes is based on finding a generalised flow g G Q{W) such 
that QifX becomes a Lie system, i.e. there exists a Lie algebra of vector fields Vq <ZV 
such that QifX G 5'(Vo). For instance, we can try to transform a particular instance 
of the systems fl24l) into a first-order differential equation associated with a nonlinear 
oscillator with a zero time-dependent angular frequency, for example, into the first-order 
system 

(25) 

= /(t)cox'". 




related to the nonlinear oscillator 



d X I 



with dr/dt = fit). 

The conditions ensuring such a transformation are 



^{t)h{t) - 7(t) = , c{t) = Co7-("+'n^), (26) 



15 



with f(t) = 7j"^(t), where 71 is a non- vanishing particular solution for 'y{t)b{t) — 7(t) = 0. 
We must emphasise that just particular solutions with 71 (0) = 1 and 71 (0) = are related 
to generalised flows in Q{W). Nevertheless, any other particular solution can be used 
also to transform a nonlinear oscillator into a Lie system as we stated. The Lie system 
( 125!) is the system associated with the time-dependent vector field 

As a consequence of the standard methods developed for the theory of Lie systems 
[13], we join two copies of the above system in order to get the first-integrals 

/. = L'2 ^ /n+i i = i 2 

h n + 1 ^ ' 

and 

^3 = ^^Hyp — — , -, 1 + 



y/h \n+r2' n+1' h{n + l) 



n + VT n + V hin + l) 

where Hyp(a, 6, c, ci) denotes the corresponding hypergeometric functions. In terms of 
the initial variables these first-integrals for g^^X read 

h = \{liit)x. - ^,{t)x,f - ^ = 1, 2, (27) 

2 7^+ {t)[n + 1) 

and 



n + 1' 2' n + 1' 7i"+i(t)/2(n + l) 



As a particular application of conditions (!26|) . we can consider the following example 
of [2S], where the time-dependent Hamiltonian 

with 7i being such that 7i(t) + uj"^ {t)'~^i(t) = 0, is studied. The Hamilton equations for 
the latter Hamiltonian are 

X = p, 

p = -sc'j;^'^'\t)x'-' - u;\t)x, ^^^^ 
being associated with the second-order differential equation for the variable x given by 

X = -sc^^;^'^^\t)x'-^ - uj\t)x. (30) 
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The latter differential equations are particular cases of our Emden equations with 

b{t) = -uj\t), c(t) = -scV^'^'H^) , n = s-l. (31) 

Notice that here the variable p plays the role of v in our theoretical development. It can 
be easily verified that these coefficients satisfy the conditions fl26l) . Therefore, we get 
that the time-dependent frequency nonlinear oscillator fl30l) can be transformed into a 
new one with zero frequency, i.e. 

^2 / 

U/ X n /„_1 

with 

T = 



dt 

wr 

reproducing the result given by Perelomov [28]. The choice of the time-dependent fre- 
quencies is such that it is possible to transform the initial time-dependent nonlinear 
oscillator into the final autonomous nonlinear oscillator. Then, we recover here such 
frequencies as a result of an integrability condition. Moreover, in view of the expressions 
( I27I) . (128!) and (13T]) . we get a, as far as we know, new time-dependent constants of the 
motion for these nonlinear oscillators. 



5.3 The Emden equation 

In this Section we apply quasi-Lie schemes to Emden equations. These equations 
appear broadly in the literature and have many applications, indeed, the review by Wong 
in 1977 [l9j contains more than one hundred cites. Here they are analysed in order to 
recover, under some integrability conditions, a set of time-dependent constants of the 
motion. 

The Emden equations we investigate are 

X = a{t)x + b{t)x'^, n=^l. (32) 

The case n = 1 is removed because it can be treated directly by means of the theory of 
Lie systems. 

Emden equations are associated with the system of first-order differential equations 

X = V, 

V = a{t)v + h{t)x''. ^^^^ 

As in the preceding examples, let us verify that the scheme S{W, V) can be used to 
handle this system. The system (1331) describes the integral curves for the time-dependent 
vector field given by 

Xt = v— + {a{t)v + h{t)x^) — , 
ox ov 

which, in terms of the basis (fT3l) for V , reads 

Xt = a(t)X4 + X3 + 6(^)^2, 



17 



so that X G S{W,V). We must remark that, as [X3,X2] ^ V, there is no Lie algebra 
of vector fields containing the vector space V" spanned by X2,X3,X4 and Xt cannot 
considered as a Lie system. 

The time-dependent change of variables induced by a control g G Q{W) transforms 
the system (155]) into 

dt \j{t) lit))'' ^ i{tf' 

dt ~ V^'^ 7(t) ait) ) ^ [ ^^a{t) a(t)7(t) a{t) ^ a{th{t) J ^ ' 

bitym 

a{t) ■ 

According to Theorem [H the latter systems describe integral curves of the time- 
dependent vector field g^^X G S{W, V). Now, we must look for a Lie algebra of vector 
fields Vq gV and a control g G Q{W) such that g^X G 5'(Vo). 

For the sake of simphcity, let us suppose that P(t) = 0. Thus, the system (IMIl leads 

to 

7(t) , a{t) , 
-X H —V , 




7(t) 7(t) 



We impose some conditions ensuring that this differential equation is a solvable Lie 
system. For instance, we want it to be of the form 



^ = fit) (Ciix' + Ci2f ') , 

dv' 

^ = fit) (C21X'" + C22V') , 

where the coefficients Cij are constant. Therefore, we get 

fit)cu = --, /(t)ci2 = — T, 

7(t) 7(t) 

/(t)c2i = fit)c22 = ait)-^ 

a[t) a[t) 



(35) 



that implies 



= 7(t) =^ -TTTT + a{t) 



cii 7W cii7(t)' 

If we fix C22 = l,cii = —1, C12 = 4 and C21 = —1, and we define Ait) = J ait)dt, the 
values of 7 and a are 



7(t) = W2 / expiAit))dt, ait) = -Lexp(A(t)) 
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with appropriately chosen constants of integration to get 7(0) = a(0) = 1 and 
g{a{t), 0, 7(t)) G Q{W). Nevertheless, any other particular solutions with different initial 
conditions can be used also. Now, since 

6(t)7"(t) la{t) 



a{t) 4 7(t)' 

we see that 

-6(t)7"+-^(t) =4exp(2A(t)) 

and 



6"^(t)exp (^^^^^-2^ j exp{A{t))dt = 0, 



which is equivalent to the expression found in [27] if we do not fix the initial conditions 
for jit) and a(t). 

Now, let us obtain a constant of the motion for ( l35l) . As a{t)/j(t) = 4/(t), we get 

exp(A(t)) 



2 / exp {A{t))dt' 

Hence, the system (l35l) admits a constant of the motion in the form 



I = -2v'^-- + x'v'. 

n + 1 



If we invert the initial change of variables, we reach the following constant of the motion 
for our initial differential equation, 

/' = |^^2 _ ^M^3,n+i^ i-2A{t)) j exp {A{t)) dt - XV exp {-A{t)) , 
which is equivalent to the one found by Sarlet and Bahar in [29] . 

5.4 Dissipative Mathews— Lakshmanan oscillators 

In this Section we provide a simple apphcation of the theory of quasi-Lie schemes to 
investigate the time-dependent dissipative Mathews-Lakshmanan oscillator 

(1 + Ax^)x - + Ax^)x - (Ax)x^ + c<;(t)x = 0, A > 0. (36) 

More specifically, we supply some integrability conditions to relate it to the Mathews- 
Lakshmanan oscillator l311 [321 [33l [38] 



(1 + Ax^)x - (Ax)x^ + A;x = 0, A > 0, (37) 

and by means of such a relation we get a, as far as we know, new time-dependent constant 
of the motion. 

Consider the system of first-order differential equation related to equation (15^ in the 
usual way, i.e. 

' i; = f , 
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and determining the integral curves for the time-dependent vector field 

Let us provide a scheme to handle the system (155]) . Consider the vector space V spanned 
by the vector fields 

X — V + ^^^"^ ^ X — ^ ^ X — fSQ) 

and the linear space W = {X3). The commutator relations 

imply that the linear spaces W, V made up a quasi-Lie scheme S(W, V). As the time- 
dependent vector field Xt reads in terms of the basis ( l39l) 



Xt = Fit)X,-uj{t)X2 + X,, 



we get that Xt E S{W,V). 

The integration of X3 shows that 

giW)^\gia{t)) = 



V — a{t)v' . 



a{t) > 0, a(0) = 1 



and the time- dependent changes of variables related to the controls of Q{W) transform 
the system (l38l) into 

x' = a(t)v', 

F(t) - ^ U' - + a{t)- 



a{t)J a(t)l + Aa;'2 ^'l + Xx'^ 
Suppose that we fix d — F{t)a = 0. Hence, the latter becomes 

x' = a{t)v', 

.. uj(t) x' , , Xx'v'"^ 

V = TT + "W- 



a(t)l+Ax'2 ''l + Ax'2' 

Let us try to search conditions for ensuring the above system to determine the integral 
curves for a time-dependent vector field of the form X{t,x) = f{t)X{x) with X E V, 
e.g. 

x' = f{t)v', 
^ =fit) T-^ + 



,1 + Ax'2 l + Ax'2 

In such a case, a{t) = f{t), Lj(t) = —a'^{t) and therefore u(t) = — exp (2 J F{t)dt). The 
time-reparametrization dr = f{t)dt transforms the previous system into the autonomous 
one 




V 



x' Xx'v 
+ 



dr 1 + Ax'2 1 + Ax'2 ■ 
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determining the integral curves for the vector field X = Xi + X2 and related to a 
Mathews-Lakshmanan oscillator (ITTI) with k = 1. The method of characteristics shows 
after a brief calculation that this system has a first-integral 



that reads in terms of the initial variables and the time as a time- dependent constant of 
the motion 



for the time-dependent dissipative Mathews-Lakshmanan oscillator fl5Bl) getting a, as far 
as we know, new t-dependent constant of the motion. 

6 Conclusions and Outlook 

We develop the theory of quasi-Lie schemes as a generalization of the theory of Lie 
systems and we prove some of their fundamental properties and find applications. In 
particular, we recover a time-dependent superposition rule for a family of dissipative 
Milne-Pinney equations. This result, which can be found in [33], is seen here from a new 
perspective as a result of systematic treatment of the family of dissipative Milne-Pinney 
equations admitting such a time- dependent superposition rule. 

Additionally, we explain from a geometric point of view some transformation prop- 
erties of time-dependent nonlinear oscillators. More precisely, we provide a geometri- 
cal understanding for some results of the Perelomov's paper [28]. Moreover, quasi-Lie 
approach allows us to investigate quantum nonlinear Hamiltonians and supply an ex- 
planation of the transformation properties for the quantum analogue of this physical 
model [34j- Finally, we also recover time-dependent constants of the motion for Emden 
equations and certain new dissipative time- dependent Mathews-Lakshmanan oscillator. 

We hope that this shows that the theory of quasi-Lie schemes can be viewed as a 
good approach to study many interesting systems of differential equations and quantum 
Hamiltonians from the same geometric viewpoint. We follow this idea in some works 
under development |35l |37] . 

Finally, the here developed examples prove that the set of time- dependent changes of 
variables allowing us to transform a differential equation in a scheme into a Lie system 
can be broader than that of those detailed for the group of such a scheme. We have 
already found an explanation for this fact which will be included in next works within 
the theory of quasi-Lie schemes [35] . 
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